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ABSTRACT 


A data acquisition and reduction system has been created 
for aircraft materials fatigue testing. The system uses an 
HP-9835 Desktop Calculator, an HP-2240A Measurement and 
Control Processor and a Materials Testing System loading 
machine. Three different computer programs on the HP-9835 
are used to analyze material properties, simulate inflight 
fatigue loading and compute fatigue damage at a stress con- 
centration. The flight loads are selected from Mil Spec 8866 
and applied in a random order. The fatigue damage at a stress 
concentration is calculated from the applied local stresses 


using Miner's Law. 
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LN ReDUeLLON 


Meee, sos 21igGh Cost Of aircraft dictates that the lifetimes 
of existing and new aircraft be extended as long as possible. 
The data acquisition system and computer programs described in 
this thesis will permit more realistic fatigue testing of 
aircraft materials; thereby enabling the engineer to predict 
and design more accurately for fatigue life. 

The programs are set up to start with a standard tension 
test to determine the modulus of elasticity and the Ramberg- 
Osgood coefficients to model the plastic region. The second 
program is for flight load simulation of a specific number 
of flight hours to be applied to a specimen with a given 
Stress concentration. This program was designed to apply 
loads that would be experienced by fighter aircraft, the 
loads being randomly ordered after being selected from the 
Military Specification 8866 1000 flight hour spectrum. The 
load data from this test is stored in sequence on magnetic 
tape, and local strain information is gathered at the point 
of stress concentration via pictures of photo-elastic 
patterns taken during the test. The last program assesses 
the damage incurred by the specimen during the flight 
Simulation. This program is specialized to predict damage 


for specimens fabricated from 7075-T6 Aluminum. 





i CONCE Es 


a. sSTRES> AND STRAIN CALCULATIONS 

Lindsey has developed [Ref. 1] a series of equations 
for calculating local stresses at stress concentration 
sites, which use a nonlinear stress-strain law and account 
for residual stresses and strains. The analysis assumes 
that the notch tip behaves in many respects like a uni- 


axial tension specimen with a Ramberg-Osgood description. 


- o gyn 
Loading the specimen of Figure 1 to a nominal stress of 


S psi, the local stress can be found by solving 


_ a aah 
KS = o{l + 8 (=) ] (2) 
where 
Ky = Elastic Stress Concentration Factor 


During the unloading portion of the cycle, Lindsey showed, 
via energy considerations, that the residual stress was 


given by 
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where o is the local stress at the notch at the peak of the 
lieadingmemet creas found from Equation (2) above. This stress 


will be compressive and will be non-zero only if os exceeds 
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eae Dreportional limit by a significant amount. The 


residual strain is positive and given by 


as! Pa oie! 
En = _= “A 8 (5) (4) 


These values of stress and strain become the starting points 
for the next cycle, and the R fatigue parameter is calculated 
taking the residual stress into account. (The R value is 

eae LacilO Of the minimum stress to the maximum stress.) The 


elastic stress in subsequent cycles is given by 


J = 0, + KS - Ee, [1 + 2 (5) 
Last 


Yield 


R 


At the next yield, a new residual stress and strain must be 


Calculated from Equations (3) and (4). 


peer ree rARY SPECIFECATIONS FOR FATIGUE 

In Order to analyze materials and structures for fatigue, 
it is necessary to have some Knowledge of the types of loads 
and the number of loading cycles that the structure will be 
Smo jected £O. For this purpose Military Specification 8866 
Flight Loading Spectrum A has been chosen (Table I). This 
spectrum was chosen by the military for loadings that were 
expected on fighter-type aircraft for each 1000 hours of 
flight time. The military expects these aircraft to have 


a lifetime of 4000-6000 hours under this type of loading. 
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TABLE [I 


PiteGHt LOAD SPECTRUM A 


3 Limit Load Number of Cycles 
35 Owe 
45 SOKO 
SS 6 ; 5010 
aS) 4,500 
75 | Be 5100 
85 Jey S50. 
S15 440 

105 LS 
die es 40 
5 16 


ee SD LeaeARD USED FOR MATERIAL PROPERTIES 

The reference material specifications for fatigue life 
were taken from NACA TN 2324 [Ref. 2]. Testing was done on 
sheet 75S-T6 Aluminum (Ai 7075-T6) for various R values. 
The results of these tests were plotted in Figure 2. 
Regression equations were formulated using the BMDP 
Biomedical Computer Programs [Ref. 3] by inputing the data 
points in Table II that were selected off Figure 2. 

The best regression formula (Table III) was found by 
using the actual R value and the stress divided by 10,000. 


The output value was the logarithm of the number of cycles 
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S-N Diagram for 7075-T6 Aluminum 


Figure 2. 
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TABLE III 


REGRESSION COEFFICIENTS 


Tea Coefficient 
INTERCEPT 1D 6452 
S moe Gan 
a 2.81098E-3 
He SSL Ola) use 
RA ~12.8099 
R® 212.476 
R® =a, Gasy0 
SR 3.68800 
s?R -0.112720 
s°R Ty eee: 
Slne 3 S865 mec 
Ee ~3.50885E-2 
stp® meGiee7E—2 
SR? Beh ao: 
S = Stress/10000 
R = Minimum Stress/Maximum Stress 
Pemrcivune, (N). Thus, the logarithm of N is the sum of the 


terms shown in Table III. This gave a multiple correlation 


coefficient of 0.96167. Table IV shows the correlation of 


the regression output to the input values of Table II. 
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Although the regression predicted output does not look like 
a close correlation to the observed data, the error of the 
regression equation falls within the scatter band for the 


Sctginal data {Ref. 2]. 
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Pete DesCREPTPEeN OF AG@UISTTION SYSTEM 


The data acquisition system consists of an HP-9835 
Desktop Computer (Controller), an HP-2240A Measurement and 
Control Processor (Processor), and a Materials Testing 
Machine (MTS). The Controller (Fig. 3) is used to generate 
the random loads, send the loads to the MTS through the 
Processor and manipulate the strain data that it receives 
back through the Processor. The Processor (Fig. 4) receives 
Gmgttal commands from the Controller to output specified 
analog values or to read analog values and send them in 
emeitaletorm to the Controller. The MTS (Fig. 5) is a 
hydraulic testing machine that applies loads specified by 
voltage Signals to a test specimen. 

The resident language in the Controller is BASIC in 
which all of the programs are written. Programs and data are 
stored on a built-in magnetic tape mass storage device. The 
Controller is linked to the Processor (or any other peripherals) 
through the HP-IB Interface Bus. This Bus has its own Read 
Only Memory ‘{ROM) and address associated with the ROM. Each 
peripheral on the Bus also has its own HP-IB address. The 
Processor is addressed by the Controller through the HP-IB 
by use of the ENTER and OUTPUT commands, which include the 


device address and the data to be sent or requested. 


69 





Figure 3. 
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HP-9835 Desktop Computer 








a 


CRATERS OU 
SS SS ; Sea 
SON 


i 
seer 


Selig et Atte? ll 





Gales. 4. HP-2240A Measurement and Control Processor 
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Figure 5. 











The Processor has four function card slots. For this 
data acquisition system one D/A card and one A/D card were 
all that were necessary. The D/A card can handle up to four 
output channels, one is used to drive the MTS and another is 
used to drive a relay, which operates an automatic camera 
shutter. The A/D card can handle up to sixteen bi-polar 
channels. In this instance eight were used by the acquisition 
program for reading strain gages, and/or load cells. Commands 
are sent to the Processor in standard ASCII format, and an 
error code is sent back to the Controller. 

The MTS is a 100,000 1lb., closed loop, feedback controlled 
Materials Test System with a hydraulic actuator. The MTS can 
be used in either load, strain or displacement control modes. 
It can be used by itself with its own controller or its own 
function generator. In this case the output from the D/A 
card was input to the MTS ina load control mode. The signal 
from the Processor was output to the hydraulic system and the 
MTS checks itself with a load cell signal until the load is 


matched. 
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IV. DESCRIPTION OF PROGRAMS 


There are three major programs written for the controller 
(in BASIC) for material properties and fatigue testing. The 
first is entitled AXIAL and is used to record data from a 
strain gaged specimen in an axial tension test to determine 
the Modulus of Elasticity, Poisson's Ratio and the Ramberg- 
Osgood coefficients for the plastic region. FLTTST, the 
second program, generates random ordered loads and outputs 
these loads through the PROCESSOR to the MTS System, which 
applies the loads to the specimen. The loads and the strains 
generated by this program are stored on tape for use later in 
data reduction. The last program, DAMAGE, takes the data 
generated by FLTTST and calculates the fatigue damage of the 


specimen. 


A. PROGRAM AXIAL 

AXIAL does not control the MTS System and so can be used 
with any system possessing voltage outputs for load and 
strain to record data. AXIAL records data while the testing 
machine performs an axial tension pull on a test specimen. 
Mic Paegcam 15 Set Uc tO read the load and both the longi- 
tudinal and transverse strains at equal intervals of time 
while the testing machine is used in a constant strain rate 
mode. The program can also be run using only the longitudi- 
nal strain gage in which case Poisson's ratio is not computed. 
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The loads are converted to stresses so that the program can 
work directly with the stress-strain relationships. After 
all the data has been taken, the program asks the user to 
specify the load at which the proportional limit seems to 
have been reached. This cutoff is entered so that a good 
linear fit can be made to find the Modulus of Elasticity (E). 
Atter the linear fit to obtain E, the value is output along 
with the slope intercept B and the correlation factor. 
The program then asks the user if a different fit is to be 
tried using a new proportional limit. If not, the program 
continues by computing Poisson's ratio and finds the Ramberg- 
Osgood coefficients using the Modulus of Elasticity previous- 
iy round. 

tiem ieadulus Of Elasticity is the most critical property 
PounaG Ine AxtAL not only for its own value in calculating 
stresses and strains but also because the Ramberg-Osgood 
coefficients are based on E. For this reason a good linear 
fit is obtained by a routine that evaluates the perpendicular 
distances of the data points from the estimated fit. The 
E that has been calculated is output to the user along with 
the correlation factor so that a judgment can be made as to 
whether the estimated proportional limit should be changed 
ieorcdommmer try tor a better fit. 

The equations governing the linear elastic portion of the 
stress-strain curve are: 

cat 
re) 
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where A is the modulus of elasticity (E). B and D are 


normally neglectibly small. Poisson's ratio is given by 


Se 
y= - G (8) 
The plastic portion of the stress-strain curve can be 
approximated by the following equation: 


a 


+ 82 


(1) 


Beta and n can be found by performing a logarithmic curve 
Piw@nretne data above the proportional limit, but the 
Ramberg-Osgood method gives as good a fit with fewer calcu- 


lations. The Ramberg-Osgood relations are: 


ecu 777) 
ee > ee (9) 
oie «(Cl ae 
g = (ht (10) 
Ar) 
where, 
¢ 5 = stress at the intersection of a line of slope .7E 
and the stress-strain curve 
O g5 = Similarly computed with a slope of .85E. 


imewes, LOr AXIAL are: 
Maximum Strain for the test 
Strain Rate to be used 
Number of Data Points desired (1000 max) 
Number of Strains to be recorded (1/2) 


Cross Sectional Area of specimen 
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Load Scale Factor (pounds/millivolt) 
Bein seale Factor (microstrain/millivolt) 


After the data run: 


Maximum Load of Elastic Range (Proportional Limit) 


me FE ROGRAM FLTTST 

Tiemmosasorltiea) program is FLTTST since it both 
controls the MTS and takes data. The key to FLTTST is the 
random selection of the load to be applied to the specimen. 
Using the Mil Spec 8866 Flight Loading Spectrum A for 
Fatigue, percentages were calculated for the number of times 
each load is to be applied (Table V). These percentages 
were then totaled and left in decimal form as numbers 
between 0 and l. Then, using the Controller's random 
number generator, the load to be applied is selected through 
a series of IF statements. For example, to select a load 
of 95% of limit load, a random number above 0.98462125 and 
below 0.99509594 has to be selected. If after going down 
through all of the IF statements a load is not selected, it 
would mean that the random number must be above 0.99961910, 
then 125 limit load has been selected, the MTS load is 
incremented in finite load steps that are user selected. 
When the desired load is reached, the Controller than 
Beemoseomemcten tle Processor read the data channels. FLTTST 
then stores these data points on tape for future reduction. 


Peeecadchneload tO be output is selected, a count is kept of 





that load so that the user can have an immediate output of 


moiombead GlStribution at the end of the run. 


TABLE V 


RANDOM NUMBER FOR SELECTING MANEUVER LOAD 


% Limit Load t Cycles Random Number 
35 ny 000 0.40470409 
45 9, 50:0 OPoes0S6226 
5D 67 5 0 0.78560206 
65 4,500 Ores 272761 
Ls 27, 500 0.95224492 
2 i, iligpges S18 Ores 26/2125 
75 440 0299509594 
105 150 0.99866686 
Eales 40 C5 ESAS TEE ALO) 


iowes tor FLITST are: 
Number of Hours to be simulated 
Number of Data Readings (maximum of 8) 
Load at One G (lbs) 
Limit Load (lbs) 
Mead scale of MTS (10,000; 20,000;50,000;100, 000) 
Load Rate (lbs/sec), and 
Data Reading Scale Factor (microstrain or pounds/ 


ime vOw tC) 
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There is another program named RDDATA that may sometimes 
be needed after running FLTTST. If data from more than one 
strain gage or load cell were recorded on the data file 
from FLTTST, it is necessary to first run RDDATA to separate 
each data file under its own heading before proceeding with 
program DAMAGE. This was done to conserve time between data 
points taken and to extend the life of data tapes. It also 
conforms with program DAMAGE since it can only compute 


damage on one set of data at a time. 


C. PROGRAM DAMAGE 

The purpose of this program is to take the loading that 
was saved on tape with the FLTTST program and calculate the 
fatigue damage accumulated at the stress concentration site 
during the flight hours simulated. The program reads the 
recorded far-field loads in the order that they were applied 
to the specimen, calculates local stresses at the stress 
riser and then uses Miner's Law for calculating the cumula- 
tive damage. Once the load creates a stress in the critical 
area greater than the proportional limit, local stresses are 
calculated using a nonlinear stress-strain law and calcula- 
tions for residual stresses and strains are made. 

The maneuver load and the one G load stored immediately 
after it are read in together. The maneuver load is then 
converted to nominal or far-field stresses and local stresses 


using the following relations: 


5 a5 (11) 


Zo 








= Kes (12) 


where, 
S = nominal stress 
L = applied load 
A = minimum cross sectional area of specimen 
o = local stress 


= stress concentration factor 


The local maneuver stress is checked against the proportion- 
al limit or the last largest plastic stress to see if it 
would cause plastic deformation. If plastic deformation is 
produced, then an iteration process must be performed to 


find the local stress, uSing the relation 


ynel 


Keoe= Ci (. lties = ] (2) 


The routine to do this uses a modified false position method 
[Ref. 4]. With the local stress determined, residual 


stresses and strains are then calculated, using 





Ih 
1 2ngE* (gy n+l ‘ 


Toe to ae UE foi) 
e) 
_ oR o\n 
Ep = zt 3 (5) (4) 


If no new plastic deformations are caused but the plastic 
region has been entered before, then local stresses are cal- 


culated using 


Ga= 90 


he settee Bep,ll 5 


Y) 
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where Siast is the largest nominal stress that has occurred 
Were | d 


which produced local stresses into the plastic region. 

Damage is then calculated using the current local stress 
and the ratio of the current stress to the stress at one G 
loading (R value). 


Using Miner's Law: 


D= N Glys)) 


where, N = number of cycles to failure and can be found using 
regression formulas that fit stress-fatigue data for 7075-T6 
aluminum. The damage caused by this cycle is then added to 
the damage caused by the other cycles computed one load at 
a time until all the data has been read. 
Inputs to DAMAGE are: 

Minimum Cross Sectional Area (sq. in.) 

Modulus of Elasticity (psi) 

Ramberg-Osgood Coefficients, and 


Proportional Limit (psi) 


onl 
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Several tests have been developed to evaluate both the 
programs and the program-acquisition system combined. AXIAL 
was checked by inputting stress and strain pairs taken from 
an existing Ramberg-Osgood representation of 7075-T6, and 
the results of the two were compared. FLTTST was checked 
by taking the output from the Processor (normally used to 
G@etve =ne MIS), routing the output from the output port to 
the input port and reading the output value on the data 
im7ewe pone. Thus, the input port was getting input just as 
fmeenommaliy would from a load cell. By doing this the actual 
voltages could be read directly from the output port with 
a voltage meter, and the input results could be treated as 
actual loads to check Program DAMAGE. This was checked by 
hand for both load spectrum and load accuracy. DAMAGE was 
checked by hand with sample runs through FLTTST as described 
above. Two full system tests were then run with the MTS, 
including photo-elastic pictures of the stress concentration 


area and strip chart recordings of the input and output loads. 


A. PROGRAM AXIAL 
The first test of AXIAL was performed by using the 
Ramberg-Osgood coefficients calculated by Kaiser [Ref. 5] to 


generate stress-strain pairs for manual input into the 


eZ 








program. The result should have been the same value for the 
modulus of elasticity and the Ramberg-Osgood coefficients as 
previously computed. Since a relatively small number of 
points were entered by hand (53) as compared to the number 
of points that will be taken by the program (1000), larger 
errors were expected on this test than will be expected on 


em actual. data run. The results of the test are shown below. 


Kaiser AXIAL 
Bee OL / EG E = 10.113E6 
n= 21.58 n = 23.475 
B = 1.479343 p= 1.517 7547 


This large difference (especially in 8) was due to the stress- 
strain pairs being calculated in stress steps of 500 psi. 

The program, using 1000 points, will have a difference of 
approximately 16 psi between points. This will make the 
accuracy of the program better than the capability of the 
system to read loads and strains accurately. 

The second test of AXIAL was an actual specimen tension 
test. A specimen of 7075-T6 Aluminum was fabricated from 
the same sheet stock as the notched panel specimens. This 
was inserted into the MTS machine which was operated in a 
constant strain rate mode. The specimen was loaded to 
pecaucera Strain of 14/000 microinches/inch, at a strain 
rate of 1000 microinches/inch/minute. One thousand data 
points were taken at constant intervals of strain during 


the test. The program then computed the modulus of elasticity 
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and the Ramberg-Osgood coefficients. These were compared 
to hand calculations made from the Load vs. Strain plot made 


by the MTS. The results of both are listed below. 


AXIAL Hane Callen lacion 
E = 9.504256 E = 9.500056 
n = 31.2154 | n = 33.0383 
8 = 3.1168E63 8 = 9.7360E66 


A plot of the actual data and the points computed using the 


calculated material properties are shown in Figure 6. 


B. PROGRAM FLTTST 

The first test performed with FLTTST was one to see if 
the random load spectrum followed and Mil Spec load spectrum. 
This was done by making several runs of 100 hour simulations 
and comparing the spectra (Table VI). As would be expected, 
the spectra did not match exactly due to the random way in 
which loads were chosen. The largest differences were in 
the area of the fewest loads, as also would be expected. 
A sample from a strip chart of the Processor output and the 
output from a load cell while the MTS is being driven by 
FLTTST is shown in Figure /7. 

Several sample tests were run to check the accuracy of 
FLTTST and DAMAGE. Knowing that a full system test was 
going to be run, the sample tests were run using the 
material properties of the specimen (Fig. 1) to be used in 
the full scale tests. These material parameters and properties 


are: 
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Figure 6. Axial Tension Test 
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TABLE VI 


SAMPLE LOAD SPECTRA 





% Load Peo LEC 
25) 17000 
45 9500 
55 6500 
65 4500 
yes 2 SH0)8. 
85 1360 
95 440 

105 510 
Pes 40 
ZS 16 
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Sample Random Loading 
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For these uses the limit load was set at the load that would 
cause local yield stress or at two-thirds of the ultimate 
stress, whichever was the lower. In this case two-thirds 
the ultimate stress was the limiting factor which gave a 
limit stress of 54,667 psi, produced by a load of 14,289 lbs. 
Assuming a 6.5-g aircraft limit load, this gave a l-g stress 
of 8409 psi, produced by a load of 2198 lbs. These are the 
limits asked for by the program FLTTST and all other loads 
are fractions of the limit load. The l-g load is used as 
the minimum load in the fatigue cycle. 

To check the accuracy and speed of different loading 
Step sizes, sample calculations were performed and compared 
with the computed actual load that should be applied (Table VII). 
For this check the material properties of the sample specimen 
were used. The conclusions were that a step size of 25 
Should be used since this gives good accuracy and takes 


approximately four minutes per hour of flight time simulated. 


C. PROGRAM DAMAGE 

Further tests were performed primarily to check the 
DAMAGE program, but the FLTTST program was used to generate 
some of the loads to run through DAMAGE. The first was a 
handmade test of ten cycles that had one load that would 
send the specimen into the plastic region locally and other 
loads that would produce recordable damage. The purpose of 


Paes Gheek was to check the program's calculation of the 
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residual stresses and strains and the difference in damage 

due to those residual stresses and strains. There were 

three loads that caused damage (85%, 95% and 125% limit load) 
and one load that caused local stresses above the propor- 
tional limit (125% limit load). Residual stresses and 

strains were calculated along with the total damage and 
compared with the program's results (Table VIII). Differences 


mcaliculations were wlthin roundoff error. 


EABEE Virtt 


TEN CYCLE LOAD SPECTRUM 





% LIMIT LOAD LOAD (1bs) STRESS R D 
5 S001 i Lotss 0.4395 0 
65 9387.08 S55 oe Ono 67) 0 
45 6430.05 24599 0.3418 0 
85 2a 5270 46465 Ore Sako eee 256 
125 iy cio tes 0 66944 0.1048 S09 SE—5 
5 5)15)(8)a es yeas, oe Sey 0 
95 5 7412-60 50488 0) Fas Ssh, 3.8419E=6 
45 6430.05 Ze 5 Oms025 0 
75 OCs SS)S\S) 7 Omaly fz 0 
55 T8565 28646 0.2448 0 
HAND CALCULATION 
ee igs psi oe 7.3442E-6 in/in 
ao oo) 
PROGRAM CALCULATION 
J, = -1313 psi lp 7.3443E-6 in/in 
De=eon oo) 7a) = 66,944 psi 


eee rrr ee 


40 





The last set of tests ran the entire system. Program 
FLTTST was used to drive the MTS, and output from the MTS 
load cell was read by the program and stored on tape. 
Photoelastic pictures were taken of the specimen before the 
test and during the test at every load above 75% of the 
mec lead. A strip chart reading of the program output 
and the load cell output (Fig. 7) were read out simultan- 
eously to enable immediate analysis of the program load 


ama MTS load correlation. 


TABLE IX 


Lone OUR ster GHi SiMULATION TEST 


+LOADS LOADS SABRC OSES: R D 
2 85 46482 Oe 77 OC Ik Oly AS pels 
5 5 SeiG2g Chalke saak 4.4607E-6 
2 105 5 ae Ora 29 IG 2B s eS 


HAND CALCULATION RESULTS 
Dea) 5.99725 —5 
PROGRAM RESULTS 


Diao. 2025 





The first test was a three flight hour simulation and 
the second was a two-hour simulation. During both tests 
the highest load applied to the specimen was 105% of the 


limit load. By reference to Table VII, 105% of the limit 
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load does not cause local stresses above the proportional 
limit; therefore, no residual stresses need to be calculated 
and all the l-g stresses will be the same. The results of 
the tests (Tables IX and X) show that the program and hand 


calculations agree within round-off error. 


TABLE X 


TWO HOUR FLIGHT SIMULATION TEST 


#LOAD SLOAD STRESS R i) 
3 85 46482 0.1770 1.0726E-6 
i 105 57576 0.1429 1.2222E-5 


HAND CALCULATION RESULTS 
D= 1.5440E-5 
PROGRAM RESULTS 


Dee 2 5234E-5 
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VI. CONCLUSIONS 


Pelomoolleved that this system will greatly simplify 
PucGmim@eredase the acctiracy of predicting fatigue damage. 
Almost all of the calculations that would normally have to 
be done by hand, one step at a time, are now performed by the 
programs. A few of the input values were left as hand 
calculations as a safety feature, so that the user would be 
aware of the magnitude of the loads at which the MTS would 
be operating. 

Although several tests and examples of program verifi- 
cation were included in this thesis, many more checks were 
performed to insure that the system was operating correctly. 
The only area in which problems occurred was with the 
Program AXIAL. If enough data points are not taken, E will 
Still be evaluated correctly but the Ramberg-Osgood coeffi- 
erenesewel! sutfer greatly in accuracy. 

Even though these programs will greatly enhance analysis 
in the fatigue testing area, the potential of the data 
acquisition system has barely been tapped. The system as 
it now stands could handle over twice as much data input 
and output. With added control cards, the system could be 


expanded at least tenfold. 
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APPENDIX A 
PROGRAM FLOW CHARTS 


PROGRAM AXIAL 
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PROGRAM FLTTST 


INPUT CONSTANTS 
AND 







CREATE DATA FILE 





> 


OQUTLOAD=OUTLOAD=- 
Sree 
SELECT RANDOM 
NUMBER SEED 


SELECT NEW 


RANDOM NUMBER OUTLOAD=lg LOAD 
















SELECT LOAD 
faa ORRESPONDS 
TO RANDOM NUMBER 


YES 
TAKE DATA AND 
STORE ON TAPE 


= 






OUTLOAD=OUTLOAD+STER 


OUTLOAD? SELECTED 
OUTPUT NUMBER 


OF EACH % LOAD 
SE ERECTED 





TAKE DATA AND 
STORE ON TAPE 
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PROGRAM DAMAGE 


INPUT 
CONSTANTS 


V 


READ MANEUVER 
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eo 

PRINT 
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